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Elasto-plastic  finite  element  analysis  has  been  used  to  model  a  single  unsized  Kevlar  49  fibre  embedded  in  an 
epoxy  matrix  subjected  to  tensile  loading.  The  predicted  axial  stress  distribution  along  the  fibre  is  compared 
with  experimental  data,  obtained  using  the  technique  of  laser  Raman  spectroscopy,  for  a  number  of 
incremental  applied  strain  levels.  The  geometry  of  the  fibre  end  has  been  modelled  employing  an  analytical 
fibre  end  fibrillation  model  (FEFM).  Furthermore,  the  influence  of  residual  stresses  in  the  system  due  to 
curing,  has  been  analysed.  Good  correlation  is  obtained  with  the  experimental  data  and  it  is  shown  that  the 
mode  of  interfacial  failure  in  the  composite  involves  yielding  at  the  fibre-matrix  interface  initiating  at  the 
fibre  end. 
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INTRODUCTION 

When  a  fibre  reinforced  composite  is  loaded,  most  of  the 
load  is  transferred  via  the  matrix  to  the  fibre.  Due  to  the 
difference  in  the  elastic  moduli  of  the  fibre  and  matrix, 
interfacial  shear  stresses  (ISS)  are  induced  at  the  fibre- 
matrix  interface.  Theoretical  predictions  of  the  ISS 
distribution  along  the  interface  have  been  made  by 
several  authors,  notably  Cox1,  for  the  case  where  both 
matrix  and  fibre  are  assumed  to  be  linear  elastic  and  the 
fibre-matrix  interface  remains  intact.  In  reality,  how¬ 
ever,  as  the  tensile  loading  on  the  specimen  increases,  the 
shear  stress  at  the  interface  reaches  a  critical  value, 
termed  the  interfacial  shear  strength  (IFSS)  for  the 
system.  Beyond  this  point  one  or  more  modes  of 
interfacial  failure  are  observed  experimentally2-9: 

•  Mode  a.  Matrix  yielding  at  the  intact  interface.  This 
mode  has  already  been  observed  to  occur  in  Kevlar 

A 

fibre-epoxy  matrix  systems  . 

•  Mode  0.  A  crack  initiates  at  the  fibre  end  and 


*  To  whom  correspondence  should  be  addressed 


propagates  along  the  interface.  Frictional  stress 
transfer  between  the  fibre  and  matrix  may  take  place 
across  the  crack  faces.  This  failure  mode  has  been 
observed  to  occur  in  glass  fibre  systems10,11. 

•  Mode  7.  A  conical  crack  initiates  at  the  fibre  end  and 
propagates  into  the  matrix  at  an  angle  to  the  axis  of 
the  fibre10'12. 

The  axial  stress  distribution  along  the  fibre  has  been 
measured  experimentally  for  crystalline  fibre-matrix 
systems  using  the  technique  of  laser  Raman  spectroscopy 
(LRS)2  9.  This  technique  relies  on  the  fact  that  the 
Raman  frequency  of  the  atomic  vibrations  of  crystalline 
reinforcing  fibres,  such  as  aramid  or  carbon,  depends  on 
the  axial  component  of  stress,  crf,  in  the  fibre  .  Unique 
Raman  frequency  versus  <jf  calibration  curves,  produced 
for  each  type  of  fibre,  are  used  to  convert  the  measured 
Raman  frequency  variation  along  an  embedded  fibre  to 
an  axial  stress  distribution. 

For  the  unsized  high  modulus  carbon  fibre-epoxy 
matrix  system,  the  axial  stress  distribution  becomes  S- 
shaped  towards  the  fibre  ends  as  the  applied  strain  is 
increased9.  The  ISS  is  zero  at  the  fibre  ends  and  has  a 
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take-up  region  followed  by  the  normal  elastic  profile.  It 
has  been  found  in  a  previous  work,  that  Mode  7 
interfacial  failure  is  consistent  with  the  experimental 
data,  for  this  carbon-epoxy  composite  system13,14. 

For  discontinuous  Kevlar  49  fibre-epoxy  com¬ 
posites2,4,6,  the  axial  stress  distribution  after  interfacial 
failure  is  trapezoidal  in  shape.  The  maximum  ISS  is 
constant  at  a  value  of  approximately  45  MPa  near  the 
fibre  ends  and  then  decays  to  zero  away  from  the  fibre 
ends,  in  a  Cox  type  elastic  manner.  The  value  of 
45  MPa  is  the  shear  yield  stress  of  the  bulk  matrix,  so 
support  is  given  to  the  micrographic  evidence  which 
indicates  that  interfacial  plastic  deformation  occurs  as 
the  applied  strain  on  the  composite  is  increased2.  This 
suggests  a  Mode  a  interfacial  failure  for  the  Kevlar 
49-epoxy  composite  system. 

Guild  et  al 4  have  modelled  a  Kevlar  49  short  fibre 
embedded  in  an  elasto-plastic  epoxy  matrix,  using  finite 
element  (FE)  analysis.  Fair  correlation  between  FE  and 
experimental  fibre  stress  data  was  achieved,  but 
exceptionally  high  von  Mises  interfacial  strains  were 
predicted  at  the  fibre  end.  These  von  Mises  matrix 
strains  were  found  to  be  much  higher  than  is  physically 
reasonable  in  a  real  composite  specimen.  For  example, 
for  an  applied  strain  of  1 .2%,  a  von  Mises  strain  of  60% 
was  predicted.  This  was  found  to  occur  for  both  a 
rounded  fibre  end  and  a  sharp  corner,  indicating  that 
the  sharp  corner  was  not  the  cause  of  the  von  Mises 
strain  concentration. 

The  reason  for  such  a  strain  concentration  was  that 
account  had  not  been  taken  of  the  real  nature  of  the  fibre 
end  geometry,  as  shown  in  Figure  1  of  Guild  et  al.4. 


When  a  Kevlar  fibre  is  cut  using  ceramic  scissors,  the  end 
is  not  cleanly  cut,  as  is  the  case  with  a  carbon  fibre,  but  is 
splayed  out  ( Figure  1).  The  cutting  process  itself  applies 
greatest  shear  stresses  to  the  surface  layers  of  the  fibre, 
which  in  the  case  of  Kevlar,  causes  extensive  fibrillation 
fo  the  fibre  end.  This  fibrillated  zone  consists  of  numerous 
Kevlar  fibrils  interspersed  with  matrix.  In  this  zone,  the 
effective  radius  of  the  fibre  is  increased  whilst  there  is  a 
corresponding  dilution  of  the  material  properties  of  the 
fibre  by  those  of  the  matrix  {Figure  2). 

To  account  for  this  complex  fibre-matrix  interaction, 
two  analytical  models  are  presented  in  which  the 
material  properties  are  inhomogeneous.  They  are  a 
one-dimensional  model,  in  which  the  material  properties 
of  the  fibrillated  zone  are  considered  to  vary  with  axial 
coordinate,  z,  and  a  two-dimensional  model,  where  the 
variation  is  both  a  function  of  axial  and  radial,  r, 
coordinates. 


FE  ANALYSIS  OF  COMPOSITES 

The  results  obtained  from  an  FE  analysis  of  a  composite 
specimen,  must  satisfy  certain  criteria.  These  criteria  are 
based  on  the  upper  bounds  on  the  values  of  stress  or 
strain,  at  any  point  in  the  specimen  model,  defined  by  the 
real  material  behaviour: 

I — At  any  point  in  the  fibre,  the  von  Mises  equivalent 
stress  must  be  lower  than  the  fracture  stress  of  the 
fibre,  Ofi 
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Figure  1  A  typical  fibrillated  Kevlar  49  fibre  end  cut  by  razor/ceramic  scissors.  The  diameter  of  the  fibre  is  12  /an 
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Figure  2  Geometry  of  fibrillated  fibre  end 


II— At  any  point  in  the  matrix,  the  von  Mises  equivalent 
stress  must  be  lower  than  the  relaxed  tensile  yield 
stress  of  the  matrix,  a^eld,  for  the  given  applied 
strain  rate: 


III— At  any  point  in  the  matrix,  the  von  Mises 
equivalent  strain  must  be  lower  than  the  tensile 
fracture  strain  of  the  matrix, 

,  * 

^  Gn 

Element  geometries  must  satisfy  the  various  mesh 
element  quality  tests,  as  regards  element  corner  angles 
and  aspect  ratios15.  Care  must  be  taken  to  ensure  that  the 
mesh  is  sufficiently  refined  in  areas  of  high  stress/strain 
gradients,  i.e.  the  discretization  error  must  be  minimal. 
This  may  be  achieved  through  the  use  of  mesh  refine¬ 
ment15,  especially  in  areas  of  high  stress  or  strain  gradients. 

The  values  of  the  various  criteria  limits,  for  the 
Kevlar-epoxy  system,  are  o\  —  3.6  GPa,  o£eld  =  31  MPa 
and  =  4 72,  for  the  given  strain  rate.  In  a  number  of 
cases16  23,  criteria  II  and  III  are  violated.  This  is 
invariably  due  to  the  fact  that  the  matrix  is  assumed  to 
behave  in  a  linear-elastic  manner.  In  addition,  some  of 
the  meshes  fail  the  FE  based  criteria,  thus  predicting 
erroneous  stress  states16,19. 


EXPERIMENTAL  MEASUREMENTS 

The  model  composite  tested  comprised  a  single  Kevlar 
fibre  embedded  in  an  epoxy  resin.  The  fibres,  supplied  by 
Du  Pont,  were  approximately  12^m  in  diameter  and  a 
Ciba-Geigy  LY5052/LY1927  epoxy  resin  was  used.  Full 
details  of  the  manufacturing  procedure  are  given  in 
ref.  4. 

The  composite  was  subjected  to  incremental  tensile 
loading.  Using  the  technique  of  LRS,  the  axial  stress 
distribution  along  the  length  of  the  fibre  was  obtained  for 
various  values  of  the  applied  strain,  eapp.  Values  of  axial 


stress  in  the  fibre  were  taken  for  each  point  along  the 
fibre,  by  averaging  the  experimental  data  at  both  ends  of 
the  fibre. 

Raman  spectra  were  obtained  using  the  514.5  nm 
(green)  line  of  an  Ar-ion  laser  as  the  excitation 
wavelength.  The  laser  beam  was  directed  through  a 
series  of  mirrors  to  a  modified  Nikon  microscope  which 
was  used  to  focus  the  laser  beam  to  a  2  ^m  spot  on  the 
fibre  via  a  suitable  objective  lens.  The  180°  backscattered 
beam  was  collected  by  the  same  microscope  objective 
and  focused  on  the  entrance  slit  of  a  SPEX  1877  triple 
monochromator  spectrometer.  In  order  to  reduce  the 
effective  scattering  volume  to  the  region  of  the  depth  of 
focus  in  the  fibre,  the  scattered  beam  was  focused  onto  a 
200  /im  pinhole  before  it  entered  the  spectrometer. 
Finally,  the  spectrometer  dispersed  light  was  directed 
to  Wright  Instruments  CCD  (charge  coupled  device) 
detector,  employed  as  the  photon  detecting  system,  and 
the  Raman  spectrum  was  recorded  on  a  PC  for  specific 
wavelength  windows.  All  frequency  peak  values  have 
been  derived  by  applying  Lorentzian  fitting  routines  to 
the  raw  data  obtained  by  the  CCD  detector. 


FINITE  ELEMENT  ANALYSIS 

The  axisymmetric  FE  model  of  the  composite  specimen 
{Figure  3a)  covers  one  quarter  of  the  diametral  section, 
due  to  the  symmetry  of  the  problem  {Figure  3b).  The 
length  of  the  fibre  considered  was  480  /zm,  which  is 
substantially  greater  than  the  transfer  length  of  the 
system.  The  mesh  consisted  of  three-noded  triangular 
and  four-noded  quadrilateral  axisymmetric  linear  ele¬ 
ments.  The  loading  took  the  form  of  a  uniform 
displacement  imposed  on  those  nodal  points  lying  on 
the  end  face  of  the  composite.  The  FE  meshes  were 
generated  using  FEMGEN/FEMVIEW  and  the  solu¬ 
tions  were  obtained  using  ABAQUS  5.2. 

In  order  to  include  the  residual  stresses  in  the 
composite  induced  thermally  during  curing,  it  is 
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Figure  3  (a)  Composite  specimen;  (b)  FE  model  geometry 


necessary  to  understand  the  method  of  fabrication. 
Matrix  resin  is  heated  to  approximately  90°C  and  the 
correct  amount  of  curing  agent  or  hardener  is  then 
added.  The  material  is  then  poured  into  a  mould  and 
after  about  15  min,  a  fibre  is  then  laid  onto  the  surface  of 
the  now  viscous  matrix.  The  fabrication  is  then 
completed  by  adding  another  layer  of  matrix.  On  cooling 
of  the  composite  specimen,  a  gelation  occurs  in  the 
matrix,  whereupon  it  goes  from  being  a  viscous  fluid  to  a 
solid.  The  temperature  at  which  this  gelation  occurs  is 
estimated  to  be  70°C. 

To  model  this  process,  a  uniform  temperature  drop  of 
50°C  was  applied  to  all  nodes  in  the  FE  model.  The 
mechanical  loading  was  subsequently  applied  to  the 
model  as  described  below.  The  thermal  expansion 
coefficient  differential  that  exists  between  the  fibre  and 
matrix  (see  Table  1)  induced  stresses  in  both  the  fibre  and 
locally  in  the  interfacial  region  of  the  matrix. 

The  material  properties  for  both  fibre  and  matrix  are 
given  in  Table  1.  The  matrix  is  assumed  to  exhibit  perfect 
plasticity  after  a  von  Mises  equivalent  strain  of  4.2%  has 
been  reached.  The  corresponding  tensile  yield  stress  at 
this  value  of  strain  was  78  MPa.  The  matrix  was 
modelled  as  isotropic  and  the  fibre  was  assumed  to  be 
elastic  and  transversely  isotropic. 

Since  polymers  exhibit  substantial  viscoelasto-plastic 
behaviour,  the  relaxed  stress-strain  behaviour  must  also 
be  considered.  The  matrix,  over  the  duration  of  the 
experiment,  undergoes  a  certain  amount  of  stress 
relaxation,  thus  reducing  by  a  certain  factor  its  effective 
elastic  modulus.  This  factor  has  been  found  to  be  0.6,  at 
room  temperature,  and  thus,  the  final  relaxed  tensile 
yield  stress  was  calculated  to  be  31  MPa2. 

An  initial  confirmation  of  the  plain  fibre  end,  as 


Table  1.  Material  properties 


Kevlar  49  fibre 

LY5052/1927 
matrix  (relaxed) 

Axial  tensile  modulus 

1 1 5  GPa 

1.78  GPa 

Transverse  modulus 

7  GPa 

1.78  GPa 

Axial  tensile  strength 

3.6 GPa  (fracture) 

31  MPa  (yield) 

Axial  fracture  strain 

2.9% 

4.2% 

Axial  Poisson’s  ratio 

0.40 

0.36 

Transverse  Poisson’s  ratio 

0.33 

0.36 

Axial  CTEa 

-5.2  x  1(T6  (K_l) 

58  x  10“6  (K-1) 

Transverse  CTEa 

41.4  x  10-6  (K“‘) 

58  x  10  6  (KT1) 

a  Coefficient  of  thermal  expansion 


modelled  by  Guild  et  al.4,  was  made.  The  fibre  was 
modelled  with  a  rounded  end,  as  shown  in  Figure  4. 

A  four  finite-size  cylinder  approximation  was  made  of 
the  fibrillated  zone  according  to  the  1-D  FEFM  (fibre 
end  fibrillation  model)  (see  equations  (Al),  Appendix). 
Each  cylinder  was  assigned  elastic  moduli  calculated 
using  a  mean  rcyl  for  the  respective  cylinder.  A  five  by 
four  discretization  of  the  fibrillated  zone  was  then 
modelled,  according  to  the  2-D  FEFM  (see  equations 
(A2),  Appendix).  Thus  each  of  the  four  cylinders, 
used  in  the  1-D  FEFM  model,  was  subdivided  into 
five  concentric  rings.  The  assignment  of  elastic 
moduli  to  each  ring  was  then  based  on  mean 
values  of  radial  distance  of  the  respective  ring.  All 
cylinders  and  rings  were  modelled  as  transversely 
isotropic,  in  line  with  the  assumptions  of  the  FEFM. 
The  levels  of  applied  strain  that  were  applied  experi¬ 
mentally  and  used  to  match  the  FE  data  to,  were  0.5, 
0.7,  1.0  and  1.3%. 


RESULTS  AND  DISCUSSION 

Simple  fibre  end  model 

The  undeformed  mesh,  in  the  region  of  the  fibre  end,  is 
shown  in  Figure  4.  The  von  Mises  strain  concentration 
which  occurs  at  the  fibre  corner  (z  =  0),  is  200%,  for  eapp 
of  1.0%,  as  shown  by  the  .von  Mises  strain  distribution 
along  the  fibre  {Figure  5).  The  reason  for  the  substantial 
increase  in  this  value,  in  comparison  with  Guild  et  al.' s4 
result  of  30%,  for  the  same  model,  is  possibly  due  to  the 
higher  degree  of  mesh  refinement  in  this  area  of  high 
strain  gradients.  A  typical  mesh  element  size  in  this 
region  was  2  pm  in  the  Guild  model,  whereas  it  is  0.4  pm 
in  the  present  analysis.  This  result  confirms  the 
inadequacy  of  the  simple  fibre  end  model  in  predicting 
the  state  of  stress  and  strain  in  the  real  fibre  composite, 
subjected  to  tensile  loading. 

Due  to  this  high  von  Mises  strain  concentration,  a 
violation  of  criterion  III  has  occurred,  thus  invalidating 
these  FE  stress  and  strain  fields  as  being  a  solution  to 
the  interfacial  failure  problem,  for  the  Kevlar-epoxy 
system. 

1-D  FEFM 

The  undeformed  and  deformed  FE  meshes  in  the 
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Figure  4  Simple  fibre  end  model:  undeformed  mesh 


Figure  5  Simple  fibre  end  model:  interfacial  von  Mises  strain  distribution  (eapp  =  1.0%) 
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Figure  6  (a)  1-D  FEFM:  undeformed  mesh;  (b)  1-D  FEFM:  deformed  mesh  (scale  factor  20)  ( eapp  -  0.5%);  (c)  1-D  FEFM:  von  Mises  strain  contour 
map  (eapp  =  0.5%) 


region  of  the  fibre  end  are  shown  in  Figures  6a  and  b .  The  areas  of  strain  concentration  are  at  the  fibre  corner  and 

fibre  is  shown  shaded  and  all  the  nodal  displacements  in  the  region  of  matrix  adjacent  to  the  fibre  end  face, 

have  been  magnified  by  a  scale  factor  of  20.  A  contour  map  The  importance  of  the  thermal  pre-load  in  reducing 

of  von  Mises  strain  (Figure  6c),  exhibits  a  maximum  at  the  von  Mises  strain  concentrations  may  be  understood  with 

fibre  comer  of  7%,  for  an  e.dpp  of  0.5%  (with  thermal  pre-  reference  to  Figure  7b,  where  the  interfacial  shear  strain 

load).  This  is  more  clearly  observed  in  the  interfacial  von  distributions  are  presented  (El 2  being  the  in-plane  shear 

Mises  strain  distributions  along  the  fibre  (Figure  7a).  The  strain).  The  thermal  pre-load  causes  a  shear  strain  of  the 
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Figure  6  Continued 


matrix  at  the  interface  ( Figure  7b),  in  a  manner  opposite 
to  the  mesh  deformations  of  Figure  6b.  When  subsequent 
applied  strains  are  imposed  on  the  thermally  prestrained 
composite,  shear  strains  opposite  to  the  thermal  shear 
strains  are  observed.  Thus,  a  shear  strain  reversal  has 
occurred.  This  effectively  means  that  the  thermal  pre¬ 
load  reduces  the  von  Mises  strain  concentration,  in  any 
subsequent  tensile  loading,  by  an  amount  approximately 
equal  to  that  generated  by  the  thermal  pre-load  alone. 

2-Z)  FEFM 

The  von  Mises  strain  contour  map  of  the  fibre  end 
after  the  application  of  thermal  pre-load  and  a  sub¬ 
sequent  applied  strain  of  0.5%  ( Figure  8),  exhibits  a 
similar  pattern  as  the  1-D  FEFM  except  that  the  von 
Mises  strain  concentration  has  reduced  in  value  to  6%. 
There  is  also  slightly  more  yielding  in  the  fibrillated  zone 
than  before.  The  stress  and  strain  fields  are  indicative  of 
Mode  a  interfacial  failure. 

The  FE  axial  stress  distributions  for  eapp  of  0.5  and 
0.7%  ( Figure  9a),  show  very  good  correlation  with  the 
experimental  data  along  the  fibre  length.  For  higher 
levels  of  eapp  of  1.0  and  1.3%  (Figure  9b),  there  is 
substantial  discrepancy  between  the  experimental  data 
and  the  FE  predicted  profiles.  The  reason  for  this 
discrepancy  is  most  probably  due  to  the  fact  that  a  more 
complex  mechanism  of  interfacial  failure  is  occurring. 
This  may  be  understood  by  considering  the  values  of 


maximum  von  Mises  strains  at  the  fibre  end.  For  the  four 
levels  of  applied  strain  of  0.5,  0.7,  1.0  and  1.3%,  the 
values  of  von  Mises  matrix  strain  maxima  are  6,  10,  37 
and  72%  (as  seen  in  Figure  7a),  respectively.  By 
substituting  these  values  into  criterion  III,  it  may  be 
seen  that  the  last  two  values  fail  the  criterion  by  a 
significant  margin,  thus  indicating  a  fracture  type  of 
failure.  The  value  of  10%  is  still  double  the  quoted 
failure  strain  of  4.2%,  but  this  is  acceptable  given  that 
thin  films  of  matrix  tend  to  fracture  at  higher  strains  than 
the  bulk  matrix.  This  is  attributed  to  the  lower  flaw 
density  in  thin  films  which  otherwise  act  as  failure- 
nucleation  sites  in  the  bulk  matrix.  A  maximum 
acceptable  value  of  matrix  failure  strain  may  be  in  the 
range  15-20%. 

It  may  thus  be  postulated  that  at  the  higher  levels  of 
applied  strain,  a  mode  of  interfacial  failure  other  than 
pure  Mode  a  is  taking  place.  Thus  it  may  be  further  said 
that  an  interfacial  failure  mode  transition  is  occurring 
from  pure  Mode  a  to  a  mixture  of  Mode  a  with  Mode  0 
or  7. 

CONCLUSIONS 

Elasto-plastic  FE  modelling  of  the  Kevlar  fibrillated  fibre 
end  geometry  has  shown  that  Mode  a  interfacial  failure 
exists  on  tensile  loading  of  a  composite  specimen,  at 
applied  strains  below  1.0%. 
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Figure  7  (a)  1-D  FEFM:  interfacial  von  Mises  strain  distributions;  (b)  1-D  FEFM:  interfacial  shear  strain  distributions 
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In  order  to  correctly  account  for  the  fibrillated  zone, 
use  of  a  fibre  end  fibrillation  model  (FEFM)  in  its  1-D  or 
2-D  forms  has  been  made.  This  model  accounts  for  the 
complex  fibre-matrix  interaction  which  occurs  in  short 
Kevlar  fibres  cut  using  ceramic  scissors. 

FE  analyses  have  shown  that  the  stress  and  strain 
fields  in  the  matrix  as  predicted  by  the  two  FEFMs  are 
similar.  The  precise  fields  may  be  important  in  trying 
to  determine  the  exact  location  of  fibre-matrix 
debonding,  which  will  inevitably  occur  at  higher 
levels  of  applied  strain,  since  the  2-D  FEFM  predicts 
a  slightly  different  stress  and  strain  field  from  the  1-D 
FEFM.  Since  little  difference  has  been  found  between 
the  stress  and  strain  fields  of  the  1-D  and  2-D  models, 
it  may  be  concluded  that  the  precise  radial  variation  of 
elastic  modulus  is  not  a  very  important  parameter. 
Hence  the  sensitivity  of  the  overall  solution  with 
respect  to  discretisation  of  the  fibrillated  zone  will  be 
minimal. 

The  residual  strains  due  to  curing  reduce  the  von 
Mises  strain  concentrations  at  the  fibre-matrix  interface 
once  the  composite  is  subsequently  loaded  in  tension. 
This  may  have  far-reaching  implications  for  the  design 
strength  of  composites,  especially  in  low  fibre  volume- 
fraction  composites  where  the  thermal  contraction  of  the 
matrix  is  significant. 

At  applied  strains  of  1.0%  or  higher,  it  is  postulated 
that  an  interfacial  failure  mode  transition  occurs  from 


pure  Mode  a  to  a  mixture  of  Mode  a  with  Mode  (3 
or  7. 

A  Mode  a  interfacial  failure  is  consistent  with  the 
experimental  data  in  the  Kevlar-epoxy  composite 
system,  at  low  applied  strains,  whereas  a  Mode  7 
interfacial  failure  has  been  found  to  be  consistent  with 
experimental  data  in  an  unsized  carbon-epoxy  sys¬ 
tem13,14.  The  reason  for  the  difference  in  interfacial 
failure  mode  has  to  do  with  the  yield  stress  of  the  matrix 
versus  the  interfacial  fracture  toughness  of  the  fibre- 
matrix  system,  on  applied  tensile  loading  of  the 
composite.  The  competing  failure  mechanisms  are 
dependent  on  fibre  geometry  and  material  property 
differentials  between  fibre  and  matrix.  Carbon  reinforced 
composites  have  a  much  higher  elastic  modulus  differ¬ 
ential  than  Kevlar  composites,  and  more  significantly, 
carbon  fibres  do  not  fibrillate  on  cutting.  This  suggests 
that  Mode  a  is  the  dominant  mode  of  interfacial  failure 
for  Kevlar  49  fibre  reinforced  composites  at  low  levels  of 
applied  strain,  whereas  carbon  fibre  reinforced  com¬ 
posites  are  more  likely  to  undergo  Mode  0  or  7 
interfacial  failure. 
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APPENDIX 

Fibre  end  fibrillation  models  (FEFM) 

To  take  account  of  the  material  property  variation  in 
the  z  direction,  we  consider  an  infinitesimally  thin 
elemental  cylinder  within  the  fibrillated  zone  {Figure  2 ), 
with  a  mean  radius  rcy]  and  thickness  dz.  The  cylinder  is 
assumed  to  contain  a  mixture  of  fibre  fully  interspersed 
with  the  matrix,  where  vf  and  vm  are  the  respective 
volume  fractions.  The  volume  of  fibre  in  the  cylinder  is 
equal  to  the  volume  of  the  fibre  in  the  unfibrillated  fibre, 
radius  rf.  Thus: 

Vf  =  {Afdz)  /  (Acy]dz) 

where  Acy]  and  Af  are  the  cross-sectional  areas  of  the 
cylinder  and  unfibrillated  fibre,  respectively.  The  fibre  is 
considered  to  be  transversely  isotropic  and  linear  elastic, 
and  the  matrix  to  be  isotropic  and  non-linear  elasto- 
plastic.  Following  a  rule  of  mixtures  argument,  the  tensile 
moduli  of  the  cylinder  are  a  weighted  average  of  those  of 
the  fibre  and  matrix,  given  by 

Ec yi(e)  =  E{  ■  v{  +  Em(e)  •  ( 1  -  v{) 


where  £cyl  is  the  cylinder  tensile  modulus  in  a  coordinate 
direction,  £f  is  the  fibre  tensile  modulus  in  a  coordinate 
direction,  Em  is  the  matrix  tensile  modulus  in  a 
coordinate  direction  and  e  is  the  given  level  of  strain. 

Substituting  for  Vf,  we  have: 

ECy\(e)  ~  Efm  {^f/^cyl}  +  ^mfr)  *  0  ”  {^f/^cyl}) 

which  simplifies  to  yield: 

£Cyl(e)  ~  £f  ‘  {rf/rcyl}~  +  £mW  *  (1  “  {rf/rcyl}2) 

Similarly,  for  the  shear  moduli,  we  have  the  correspond¬ 
ing  expression: 

Gcyci(e)  =  Gf  •  {rfAcy,}2  +  Gm(e)  •  (1  -  {rf/rcyl}2) 

where  Gcy\  is  the  cylinder  shear  modulus  in  a  coordinate 
plane,  Gf  is  the  fibre  shear  modulus  in  a  coordinate  plane 
and  Gm  is  the  matrix  shear  modulus  in  a  coordinate 
plane. 

Since  the  second  term  in  each  equation  is  always  much 
smaller  than  the  first  term,  over  the  range  of  strain 
considered,  we  may  neglect  the  non-linearity  of  the 
matrix  stress-strain  curve,  without  much  error.  This 
simplifies  the  equations  to  give: 

ECyl  ~  £f*  {rt/rcy\Y  +  Em  *  (1  ~~  {rf/rcyl}2) 

Gcyl  =  Gf  •  {r f/ rcyl}2  +  Gm  •  ( 1  -  {rf/rcyI}2)  (Al) 

Both  equations  define  the  variation  for  all  the  various 
elastic  moduli,  since  the  fibrillated  zone  is  transversely 
isotropic  due  to  the  influence  of  the  fibre.  Hence,  the 
elastic  moduli  of  a  given  cylinder  are  solely  functions  of 
axial  coordinate,  z,  along  the  fibrillated  zone;  as  z 
decreases,  rcyl  increases  linearly.  When  rf  =  rcyi  the 
equations  reduce  to  the  desired  result  for  the  fibre  elastic 
moduli. 

In  the  2-D  FEFM,  we  consider  that  the  material 
properties  vary  in  both  the  axial  and  radial  directions.  It 
is  probable  that  towards  the  axis  of  the  fibrillated  zone, 
penetration  of  the  matrix  becomes  progressively  more 
difficult.  The  extent  of  fibril  dilution  by  the  matrix  will 
therefore  be  maximum  at  the  periphery  of  any  given 
cylinder,  and  minimum  at  the  axis.  Thus,  the  material  in 
the  peripheral  region  of  the  fibrillated  zone  will  have 
lower  elastic  moduli  than  that  nearer  the  axis.  It  is 
assumed  that  the  variation  of  material  properties  is  linear 
with  radial  axis  coordinate,  and  also  that  the  elastic 
moduli  at  the  axis  of  the  zone  are  three  times  their 
corresponding  values  at  the  periphery.  The  average  value 
of  the  radially  varying  tensile  elastic  modulus  of  a  given 
elemental  cylinder,  Ecr(r),  is  equal  to  £cyl.  Hence,  we 
have 

Ect  ( r )  =  ar  +  b 

where  a  and  b  are  coefficients  determined  by  the 
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boundary  conditions: 


moduli  variations 


£cr(0)  =  1.5  Ecyl  Ec[(rcy] )  =  0.5  Ecyl 
Solving  for  a  and  6,  we  have  for  the  tensile  and  shear 


Ecr(r)  =  Ec yi  { 3/2  -  r/rcyl} 
Gcr(r)  =  Gcy ,{3/2  -  r/rcy\} 
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